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ABSTRACT 


Monitoring soil moisture changes by remote means requires timely 
and repetitive sensor coverage. This dictates that the sensor must 
have all-weather capabilities. Microwave sensors have this capabil- 
ity. Until recently, microwave remote sensing of soil moisture has 
been limited to determining the best single sensor with the maximum 
sensitivity to soil moisture. Visible/near-infrared sensors have been 
developed that estimate biomass and discriminate between crops. Vege- 
tation over a particular area does not change as rapidly as surface 
soil moisture. Therefore the absence of all-weather capabilities is 
not a severe limitation when visible/near-infrared biomass estimates 
are combined with microwave sensors to estimate soil moisture. 

MuUi frequency sensor data from NASA's C-130 aircraft were used 
1) to determine which of the all-weather microwave sensors demon- 
strated the highest correlation to surface soil moisture over optimal 
bare soil conditions, and 2) to develop and test techniques which use 
visible/infrared sensors to compensate for the vegetation effect in 
this sensor's response to soil moisture. Soil moisture and aircraft 
data were collected at two agriculture’ areas; Go, i, Oklahoma in 
1978 and Dalhart, Texas in 1980. 

The L-band passive microwave radiometer was found to oe the most 
suitable single sensor system to estimate soil moisture over bare 
fields. In comparison to other active and passive microwave sensors 
the l-band radiometer 1) w^s influenced least by ranges in surface 
roughness, 2) demonstrated the most sensitivity to soil moisture 
differences in terms of the range of return from the full range of 


soil moisture, and 3) was less sensitive to errors in measurement in 
relation to the range of sensor response to soil moisture that the 
active microwave systems. L-band emissivity related more strongly to 
soil moisture when moisture was expressed as percent of field 
capacity. The perpendicular vegetation index (PVI) as determined from 
the visible/infrared sensors was found useful as a measure of the 
vegetation effect on the L-band radiometer response to soil moisture. 
A linear equation was developed to estimate percent field capacity as 
a function of L-band emissivity and the vegetation ifdex. The predic- 
tion algorithm improves the estimation of moisture significantly over 
predictions from L-band emissivity alone. 


INTRODUCTION 


Estimation of soil moisture over large areas is useful in meteor- 
ology, hydrology, and agriculture. In meteorology, atmospheric models 
require information about the energy flux at the earth's surface. The 
two types of energy exchanged at the surface are sensible heat and 
latent heat. Sensible heat absorbed and released by the soil is for 
the most part, a small component of the surface energy balance. When 
latent heat is considered, some measurement of soil wetriess is needed 
in order to relate actual evaporation to potential evaporation rates 
(Priestley and Taylor, 1972). In hydrology, real-time estimates of 
soil moisture condition at the beginning of a storm event would 
improve the ability to estimate runoff and provide flood warnings that 
would save both lives and property. In agriculture, a time history of 
surface soil moisture can be used to determine the severity and areal 
extent of drought conditions as well as inputs into soil moisture 
profile models needed to estimate agricultural yields. This 
information is useful to monitor the food supply for the growing 
population of the world. 

Until the advent of le computer and space age technology there 
was no reasonable way to handle the volume of data required to map 
soil moisture estimates over an area as large as the United States, 
much less on a global basis. Conventional soil moisture measurement 
techniques are very time-consuming and therefore measurements are not 
widely or regularly obtained over most of the United States. The 
spatial variations of soil moisture make it difficult to extrapolate 
conventional point measurements to represent an integrated value over 


a large area. These restrictions have led to research in developing 
remote sensing techniques to estimate soil moisture. The remote sens- 
ing approach to soil moisture estimation requires a new scale of 
thought. Sensing devices, particularly when satellite mounted, inte- 
grate many of the micro-scale variations that have been noticeable in 
point measurements. However, these remotely sensed measurements may 
actually give the operational user more meaningful information than 
point measurements, when used with a new set of models which use 
remote sensor inputs (Blanchard, 1979). 

To monitor surface soil moisture by remote means, timely and 
repetitive sensor coverage is required. Soil moisture changes rapidly 
in the near-surface layers. After an irrigation or rain, the soil 
surface dries out within a few days in clear sunny weather (Taylor and 
Ashcroft, 1972). This dictates that the sensor must have all-weather 
capabilities. Visible and infrared sensors do not meet the all- 
weather criterion because the short wavelengths in this region of the 
spectrum are attenuated and many times completely masked by cloud 
cover. Microwave sensors with wavelengths greater than 4 cm are not 
hindered by cloud cover or lack of sunlight. Also, longer wavelengths 
penetrate into the soil medium to a certain extent (Lundien, 1966). 
This provides microwave sensors with two potential advantages over 
visible and infrared sensors: 1) all-weather capability, and 2) the 
ability to gather information about the subsurface as well as surface 
moisture conditions. 

At visible and infrared frequencies, the sensor response contains 
information only about a very thin layer at the air-soil interface 
(Ulaby, 1974). However, this does not mean that visible and infrared 
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sensors are not valuable tools for large area soil moisture monitor- 
ing. Techniques using visible/infrared data have been developed that 
discriminate between crop types and estimate biomass, leaf area index, 
and percent cover (Richardson and Wiegand, 1977 and Rouse ^ al . , 
1973). Jackson ^aj_. (1981b) suggested the use of vegetative biomas' 
estimates from visible/near-infrared measurements in combination with 
microwave sensors to determine soil moisture. Vegetation over a 
particular area does not change as rapidly as surface soil moisture. 
Only three or four scenes in the visible and near-infrared wavelength 
region of the spectrum may be necessary to classify the vegetation 
during a season. Therefore the absence of all-weather capabilities is 
not a severe limitation when visible/infrared biomass estimates are 
combined with microwave sensors estimates of soil moisture. 

Three basic terrain parameters affect microwave sensors: 1) soil 
moisture, 2) roughness, and 3) vegetation. The frequency and incident 
angle determine the scale of roughness to which microwave sensors are 
sensitive. The influence of vegetation is also dependent upon fre- 
quency and incident angle. Until recently, research in microwave 
remote sensing of soil moisture has been limited to determining the 
best single sensor which maximizes sensitivity to soil moisture while 
minimizing the effects of rougnness and vegetation. However, the 
optimum frequency and incident angle may change for each agricultural 
region with different tillage practices and types of crops that are 
grown. 

Active microwave investigations at the University of Kansas have 
identified the 4.25 GHz active microwave with incident angles near 10 
degrees as the one sensor that can measure soil moisture over all 
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terrain conditions (Ulaby and Batlivala» 1976; Ulaby ^ 1978). 

However, data presented by Ulaby ^ (1977) showed a 3 dB shift in 

the sensor response and a reduction in sensitivity caused by rough- 
ness. Active microwave has also shown promise in crop classification 
when used at large incident angles where the effects of roughness and 
vegetation dominate (Bush and Ulaby, 1978; Schwarz and Caspall, 1963; 
El lermeier ^ , 1967). 

Passive microwave sensors are also affected by the terrain para- 
meters. Numerous investigations, such as Newton (1977), Schnugge 
(1978, 1980a, 1980b, and 1980c), and McFarland (1976) have demon- 
strated the sensitivity of passive microwave sensors to soil mois- 
ture. Lee (1974) and Newton (1977) demonstrated that surface rough- 
ness influences this sensitivity to soil moisture. The effect of 
vegetation on passive microwave response was identified in studies 
such as Lee (1974) and Ki rdiashev ^ j^. (1979). 

Very few investigations have combined the responses from various 
sensors to estimate soil moisture. Using L-band passive microwave 
system, Newton (1977) used the response from two polarizations at one 
incident angle to infer the roughness contribution on the microwave 
response at a different incident angle. Investigations which combine 
microwave sensors with sensors from other regions cf the spectrum have 
not been performed because concurrent data have not been available. 

The attention of this study is focused on investigating techniques 
which infer soil moisture over agricultural lands from a combination 
of aircraft-mounted sensing devices. The sensors included visible, 
infrared, and microwave. Multi-sensor data over a wide rang? of 
frequencies and incident angles are available only from NASA's C-13C 
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aircraft. Two of the most comprehensive soil moisture experiments 
using this aircraft were performed jointly by Texas A&M University and 
the University of California at Santa Barbara at sites in Quymon, 
Oklahoma and Oalhart, Texas. These data were used in this study 

1) to determine which of the all-weather microwave sensors 
demonstrates the highest correlation to surface soil moisture over 
optimal bare soil conditions. 

2) to develop and test techniques which use visible/infrared sensors 
to compensate for the vegetative effect in the microwave sensors' 
response to soil moisture. 

It is hypothesized that the final technique would have the potential 
for estimating soil moisture over a variety of agricultural crops with 
only remote sensing inputs. These techniques could then be expanded 
to larger areas with a satellite mounted sensor package. 

PHYSICS OF MICROWAVE INTERACTIONS 

There is some conceptual difficulty in relating to the microwave 
portion of the spectrum. Visible and near-infrared sensors are 
extensions of our eyes. Thermal sensors can be identified with 
walking barefoot over hot pavement or cool, lush grass. However, we 
cannot put our hand on a surface and feel the parameters which 
influence microwave response. The emission, transmission, scattering, 
and absorption of microwave radiation are governed by the same 
physical laws that govern ladiation in the visible and infrared 
portions of the spectrum. Microwave energy penetrates through the 
atmosphere and clouds with minimal attenuation, giving the microwave 
sensors all-weather capabilities. Since microwave energy has longer 
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wavelengths, they are affected by parameters below the surface as well 
as those on the surface. This affords the opportunity to gather 
information about the soil volume itself. 

The ability of microwave sensors to detect soil moisture is based 
upon the difference in dielectric constant of dry soil and of water. 

t is common practice to express the dielectric constant, k, in terms 
of the dielectric constant of free space, kq* 

Kp= <0 ic (1) 

where is a dimensionless complex function expressed by 

Kp = <p' - jKp" (2) 

where <p' is the relative permittivity, is related to conduc- 

tivity, and j denotes the complex portion of the function. The magni- 
tudes of both Kp' and <p" are much greater for water than for dry 

soils and air space at microwave frequencies. In a dry soil, the 
value ot Kp* is typically between 3 and 5, whereas the value of 
ter' ’S about 80 for water (Schmugge ^ , 1974). 

Investigations such as those by Cihlar and Ulaby (1974), Hoekstra 
and Delaney (1974) and Newton (1977) have presented the relationship 
between p'^rmitti vity (tcp* end icp") and moisture content for 
several different soil-texture types and frequencies. An example of 
the texture effects at 1.4 GHz is presented in Fig. 1 (Newton, 1977). 
Both the sand and clay soils have approximately the same value of 
<r' at 0% moisture. The function of tcp' for each soil has two 

regions separated at a "transition" value of water content. For sand 
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the transition value is approximately 3% moisture content in contrast 
to 20% for the clay. Newton (1977) attributes the change of Kp' 
below the transition zone with the increase of the tightly bound or 
hygroscopic water content. As more water is added, the more loosely 
bound pore water causes a larger effect (slope). This shift agrees 
with the established principle that with increasing clay content there 
is an increasing percentage of tightly bound unavailable water (Brady, 
1978). Newton (1977) concludes that the dielectric properties are 
directly related to availability of soil water to vegetation indepen- 
dent of soil type. It is important to note that <p is a function of 
frequency. 

Microwave sensors are of two basic types, active and passive. 
Active microwave sensors (radar and scatterometers) illuminate the 
target with microwave energy and measure the amount of energy returned 
to the antenna. Passive microwave sensors (radiometers) measure the 
amount of microwave energy that the target emits. 

The terrain-sensor interaction processes in microwave remote 
sensing are governed by the target geometry and dielectric proper- 
ties. These properties depend upon their relation to the sensor 
parameters (Cihlar and Ulaby, 1974). Sensor parameters include 
frequency, incident angle, and polarization as well as whether the 
sensor is active or passive. The ability of microwave sensors to 
detect soil moisture is obtained from the direct relationshir between 
microwave response and the dielectric properties of the soil (Hoekstra 
and Delaney, 1974). As the soil moisture changes, its dielectric 
properties also change. This causes a corresponding change in the 
microwave response. 
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Target geometry refers to the slope, aspect, degree of roughness, 
and type and extent of vegetative cover (biomass). The effects of 
target parameters are all integrated into one microwave response. 
Consequently, it is difficult to relate physical measurements of 
target parameters to the effects they have on microwave responses. 
This is further complicated because these effects are also a function 
of the sensor parameters (Rouse 1969). 

Active Microwave 

Radar return is that portion of the transmitted energy which 
returns to the radar receiver. The radar cross section, o, or alter- 
nately the radar cross section per unit area, (scattering coefficient) 
a°, depends upon the properties of the target. Thisparameter contains 
all the information about the illuminated terrain that a radar is 
capable of sensing. The scattering cross section, am> a single 
point scatterer is related to the radar return, by inverting the radar 
equation and solving for om (Rouse 1969): 


a 

m 


K, (4.)V 


— “yT" 

t m 


( 3 ) 


where: 

Om = scattering cross section for a single point scatterer 
(m2) 

Pp = received signal power (watts) 

P^ * transmitted signal power (watts) 

Dm = slant range distance to the mth scatterer (m) 

A = wavelength of the transmitted signal (m) 
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* antenna gain (dimensionless) 

The average scattering coefficient per unit area for many scatterers 
Is a dimensionless real numoer whose magnitude Is a function of the 
terrain parameters. Since o° Is dimensionless It can be analyzed In 
dB and expressed as a function of terrain parameters: 

0 ° » f(x, e, P, K, r) (4) 

where 

0 = angle of Incidence relative to nadir 
^ - aspect angle 

P = polarization of the Incident wave 
K = complex dielectric constant 
r = Index of surface roughness. 

An analytical expression for is not readily obtainable without 
simplifying assumptions (Rouse 1969). 

In the following discussion, a perfectly smooth surface at the 
air-soil Interface Is assumed to facilitate the understanding of the 
relationship of o° to soil moisture. Under these conditions the 
Fresnel reflection coefficient, r, can be defined. This coefficient 
relates the magnitude and phase of the reflected electric and magnetic 
fields to those of the Incident fields of a plane wave which strikes 
the interface. The parameter of interest is the power reflection 
coefficient, R: 


R 



( 5 ) 




( 6 ) 
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OKIUtNAL ^AQS 53 

OF POOR QUALITY 


[{pK^' " P)^+ (wic.'' - q)^] 

Ry = 2 7 (7) 

^ [(u</ * pf* iuK," - qr] 

where H and V represent horizontal and vertical polarization and y * 
cose. The parameters p and q are given by: 

P lll'r' * “ - • fv' * (8) 

P ■ 7T - (V" * - ‘>1*^^ (8) 

The sensitivity of R to the complex dielectric constant is not readily 

apparent from these equations. Using approximate values of Kr 
Fig. lb and equation (6), for a range of look angles is presented 
in Fig. 2 . Small incident angles show the most sensitivity to changes 
in soil moisture. The term "look angle" is synonymous with incident 
angle. 

The Frtsnel approach addresses only smooth or specular return 
from the surface with no subsurface contribution or depolarization. 
The energy is reflected as from a mirror, with the angle of incidence 
equal to the angle of reflection. Spectral reflection is seldom pres- 
ent in natural terrain except for areas with smooth water surfaces. 

At the other extreme, isotropic scatterering is when the scattering 

coefficient is independent of the angle of incidence. The "degree of 
roughness" of the terrain, a relative quantity, dictates the extent of 
scattering which forms the reradiation pattern. Examples of specular 
and isotropic scattering along with an intermediate case are presented 
in Fig. 3. 
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LOOK ANGLE IN DEGREES 


FIG. 2 Relationship between power reflection coefficient (Rh) and 
look angle for 0% and 30X soil moisture. 
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FIG. 3 Examples of spectral and Isotropic 
Intermediate case. Incident angle 


scattering along with an 
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A large number of scattering theories havt been developed by 
Barrick (1968), Kodls (1966), Fing (1966), Beckmann (1965), Semenov 
(1966), Beckmann and Spizzichino (1963), Leader and Dalton (1972), 
Rouse (1972) and Blanchard (1977). The authors have modeled the 
effects of surface roughness, subsurface contribution, and 
depolarization. These modeling efforts have had varyi'*'’ degrees of 
success for particular terrain roughnesies but no one theory works 
universally over all scales and distributions of roughnesses and 
subsurface conditions. 

Vegetation effects further complicate the dependency of o° on 
soil moisture. Vegetation scattering models have been developed by 
Peak (1959a), Du (1969), and Waite and Cook (1974) with varying 
degrees of success. The physical phenomena of radar scattering is 
very complex and not well understood. The scope of this study does 
not allow in-depth investigations of radar reflectivity theory. Such 
a study must be done on small, controlled surface arear with accurate 
measurements of roughness, vegetation, and soil dielectrics. However, 
it important to recognize the impact of scattering from surface 
roughness and vegetation on the scattering coefficient in the attempt 
to relate the scattering coefficient to soil moisture. 

Passive Microwave 

Since the microwave emission of a soil originates in the so’l 
volume and propagates outwardly, it is reasonable to examine the emis- 
sion process within the volume separately from the modifying effects 
of surface roughness and vegetative cover. The fundamental concepts 
of microwave radiation are presented here. These concepts are des- 
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cribed in more detail by Paris (1969 and 1971), Marlon (1975), Newton 
(1977), Hess (1959), and Linti and Simonett (1976). 

The specific intensity of radiant energy emitted by a black body 
is expressed by Planck's radiation law: 

®bb == ^ (exp(hv/CT)-lf^ (10) 

X 


where 

8bb " specific intensity of radiant energy emitted in a 

particular frequency range that passes through an element 
of area in an interval of time and surrounded by a solid 
angle (watts m’^ steradian"^ Hz*!) 

C - Boltzmann's constant (1.38 x 10*^3 JK'!) 

f 1 tn '*»rature of the blackbody (K) 

X = wavelength (tti) 

h = Planck's constant (6.623 x lO*^'! Js) 

V ' frequency (Hz) 

At microwave frequencies and surface temperatures typical of 
Earth, the Ra.vleigh-Jeans approximation of Planck's e<iuation for 
thermal radiation fron; a blackbody applies. A Maclaurin expansion of 
the bracketed quantity in equation (10) is: 

exp(hv/CT)-l » hv/CT ♦ h^ v^/2t‘'/T^+ ... (11) 

Equation (ID) reduces to the Rayleigh-Jeans approximation 


B 


bb 


2CT 


( 12 ) 


IS 


when the higher order terms in equation (11) are neglected as very 
stna 1 1 . 

The emission from a blackbody is completely random; therefore, 
the polarized specific intensity of emission (Bp) can be expressed 
as 


(13) 

since one-half the total intensity of emission would exist in each 
plane of polarization. Material objects or, the earth's surface are 
not perfect radiators, so they are sonietimes referred to as gray- 
bodies. The intensity of graybody radiation may be thought of as the 
intensity radiated by an equivalent blackbody multiplied by an effici- 
ency factor. It is assumed that the form of the emission is the same 
as that of a blackbody and that t.,e polarized intensity. Ip, is 
given by 

Ip = (14) 

where t is the polarized emissivity of the gray body. 

The polarized brightness ter^>erature, Igy, is defined by 

Trt = <T> 'p 

which reduces to 

Tgy = eT (16) 

10 






where T is the teti^erature of the emitting graybody. Simply stated, 
all objects with temperatures above absolute zero emit electromagnetic 
energy due to the thermal agitation of charged particles. The radia- 
tion intensity is linearly related to the temperature and physical 
properties of the object. 

Various sources contribute to the radiation received by a passive 
microwave system. The most important contribution for this study is 
the contribution from the earth's surface (cTsurf)* Other major 
sources include the atmosphere between the antenna and the surface, 
Tatm» downward radiation from both sun and the atmosphere that 
is reflected fron; the surface to the radiometer, Tsky* For a 
constant height, frequency, and angle 

^BT ^ATM ^ ^ ^SURF^ 

where 

Tgj = brightness temperature, as seen by the radiometer (K) 

T = atmospheric transmission (dimensionless) 

1-e = reflectance of the surface (dimensionless) 
e = emissivity of the surface (dimensionless) 

T^urf thermometric temperature of the surface layer (K) 

Most studies, including this one, neglect all atmospheric effects 
(^ATM* ^SKY=0, T-1). This simplifies equation (17) to 

^BT ^ ^SURF* 


17 


This assumption is valid for the majority of atmospheric conditions, 
but as shown by Paris (1971), it loses validity in the presence of 
relatively large particles such as rain drops in the atmosphere. 

When applied to a soil this simple approach to microwave emission 
assumes a semi -inf inite soil medium with a smooth surface that is in 
thermodynamic equilibrium and exhibits blackbody (isotropic) radia- 
tion. Peake (1959b) demonstrated that e is equal to absorptivity. 
This is the basis for describing the emission from a soil volume, cq, 
in terms of the power reflection coefficient; 

eo = 1 - «p (19) 

where p denotes polarization and Rp can be calculated from equation 
(6) or equation (7). In Peake's analysis uniform soil moisture and 

temperature distributions with depth were assumed. Richerson (1971) 
and Casey (1972) used an effective Fresnel coefficient from horizontal 
homogeneous layers to provide for non-uniform moisture profiles. 
Strogryn (1970) developed a rigorous solution for non-uniform 
temperature and moisture profiles using Maxwell's equations. Using an 
approach similar to Strogryn, Tsang jl_. (1975) formulated emission 
from a stratified medium. Numerical results were illustrated and 
compared with the closed form analytical solutions. 

Radiative transfer techniques (Newton, 1977; Burke and Paris, 
1975) provided a simpler alternative to the general formulation of 
Stogryn (1970). The soil volume was modeled as a horizontally strati- 
fied media with each layer consisting of non-scattering homogeneous 
soil. Soil moisture and temperature were assumed constant in each 
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layer. This approach lends itself to physical interpretation and is 
not computationally difficult. 

In nature, the soil surface is seldom perfectly smooth and with- 
out vegetation. The remainder of this section examines the influence 
of surface roughness and vegetation upon the passive microwave res- 
ponse to soil moisture. 

Surface roughness modifies the microwave emission from a soil 
volume. Choudhury ^ £l^. (1979) modeled roughness by modifying 

equation (19) with a roughness parameter, h: 

e = cq + Ac (20) 


where 


Ac = Rp(l - e"^ (21) 

Using both 1.4 GHz and 10.6 GHz radiometers, Lee (1974) and Newton 
(1977) reported measurements of bare, smooth, and rough soil surfaces 
that demonstrated that increased surface roughness decreases the sen- 
sitivity of the emission to changes in soil moisture, and that for dry 
soils the emission is approximately independent of surface roughness. 
Using these same frequencies, Blinn and Quade (1974) obtained measure- 
ments for both raked and smooth sand. From their results, the 10.6 
GHz frequency was affected much more by that specific scale of 
roughness than the 1.4 Ghz frequency. This indicates that the 
"sensed" roughness is a function of wavelength. For a smooth field 
Newton (1977) reported a 100 K change in Tbj in going from wet to 
dry soils. This sensitivity was reduced to about 60 K for extremely 
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rough fields that had been deep plowed to produce large clods. This 
shows good agreement with the 80 K range of Igy for the moderately 
rough fields encountered in Phoenix (Schraugge, 1978). 

The effects of vegetation on the surface brightness temperature 
have been modeled by Sibley (1973), Basharinov and Shutko (1978), and 
Kirdiashev ^ aj_. (1979). Experimental results from Lee (1974) sup- 
port Sibley's model which predicts that the vegetation is essentially 
an attenuator at low vegetation density and a predominant emitter at 
high vegetation density. As an attenuator, the vegetation lowers the 
apparent temperature measurement; as a predominant emitter, it raises 
the apparent temperature. Lee (1974) used 1.4 GHz and 10.6 GHz 

frequency radiometers and reported that for 10.6 GHz, thick vegetation 
had completely masked the soil contribution regardless of surface 
roughness. In contrast, for the 1.4 GHz frequency, the vegetation 
generally exhibits an attenuating effect. Therefore, the effect of 
vegetation on apparent temperature is dependent on frequency, with 
lower frequencies showing more potential to detect soil moisture 
through the vegetation. 

REVIEW OF PREVIOUS WORK 
Active Microwave 

Early ground-based measurement programs at Ohio State University 
(Taylor, 1959 and Cosgriff ^ a_[., 1960) and by the Corps of Engineer 
(Lundien, 1966) investigated radar responses and penetration into a 
variety of natural and man-made surfaces. Roughness, vegetation, and 
soil moisture were found to be the general classes of parameters which 
influenced radar response. Due to the we 1 1 -developed photographic 
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interpretation techniques, the greatest excitenent in the remote 
sensing community during the 1960s and early 1970s was created by 
radar image analysis. Numerous authors presented qualitative image 
interpretations of the effects of vegetation (Ellermeier ^ , 1967; 

Moore and Simonett, 1967; Morain and Simonett, 1966; Schwarz and 
Caspall, 1968; and Berger, 1970) and suggested radar as a crop 
classifier. However, the launching of Landsat redirected vegetative 
studies away from radar and toward visible/infrared tevhniques. 

MacDonald and Waite (1971) renewed interest in radar's potential 
for soil moisture detection. Their study was the first to qualita- 
tively discern soil moisture variations in radar imagery at relatively 
high frequencies and large incident angles. Previous experiments with 
similar imaging radars had been conducted with soil moisture showing 
little influence upon radar return even for bare soils (Simonett, 
1970). However, MacDonald and Waite (1971) discerned the difference 
between wet swampy lands and the better drained, relatively drier 
soils of natural levees. Both these areas were beneath defoliated 
tree branch cover. Studying flooded low areas in Central Texas, 
Blanchard (1979) described this phenomena as multiple scatter- 

ing between the trees and the smooth free water surface. 

The University of Kansas has since provided the bulk of research 
into the dependency of radar response on soil moisture. A description 
of the truck -mounted Microwave Active Spectrometer (MAS) system is 
provided in Ulaby ^ aj,. (1978). The research using this system iden- 
tifies C-band radar (4.25 GHz) at incident angles near 10 degrees as 
the optimum frequency which can estimate soil moisture over all ter- 
rain conditions with a minimum of influence from roughness (Ulaby and 
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Batlivala, 1976; Ulaby^jl^., 1978). However, sensitivity (approxi- 
mately 3 dB) to roughness (Ulaby ^ £]^. 1977) and vegetation type 
(Ulaby^^,, 1979) is demonstrated at this frequency and incident 
angle. They found that as soil moisture varies from 0 to 30%, o° 
changes by approximately 10 dB. Jackson ^ (1981a) reported 

agreement with the University of Kansas with aircraft mounted radars 
by illustrating that C-band at 15 degrees incident angle exhibited the 
best linear correlation with soil moisture. However, the 14 points 
presented exhibited a sensitivity of over 20 dB per 30% moisture 
range. 

Blanchard (1978) presented results from aircraft-mounted radar 
experiments (13.3 and 1.6 GHz). The influence of roughness and 
vegetation were minimized at angles between 10 and 20 degrees. At the 
10 degree incident angle, the moisture sensitivity increased from 4 dB 
at 13.3 GHz to 9 dB at 1.6 GHz per 30% soil moisture. This closely 
agrees with the Kansas data for 4.26 GHz reported by Ulaby 
(1977) and indicates no significant difference in sensitivity between 
4.25 and 1.6 GHz. However, the 1.6 GHz data from Blanchard (1978) 
does not agree with the low sensitivity shown by Ulaby at 1.1 GHz. 

Most radar measurements are accurate to only 3 dB. As mentioned 
before, roughness has at least a 3 dS effect at the optimum incident 
angles. The magnitude of these variations along with vegetation 
effects can significantly mask the effectiveness of radar to quantify 
soil moisture when applied to a variety of terrains. 


Passive NIcroMave 


Many investigators have experimentally shown the sensitivity of 
truck and ai rcraft -mounted passive microwave radiometers to soil 
moisture variations. Poe ^ al^. (1971) observed that the emissivity 
of a smooth bare field varied fran 0.5 for very wet soil to greater 
than 0.9 for dry soil. Blanchard (1972) and Blinn and Quade (1972) 
reported decreasing Igy with increasing moisture content. The rates 
ranged from 1.5 to 4.8 K per one percent change in moisture content 
depending on frequency and soil type. This range is supported by 
Basharinov et (1974), who reported a slope of approximately 3-4 K 
per one percent moisture change for the 1-3 cm wavelength range. A 
comparison of aircraft observation of Tpy over agricultural fields 
around Phoenix in 1973 (Schmugge ^ aj^. , 1974) and 1975 (Schmugge, 
1976; Choudhury , 1979) demonstrated that these results are 

repeatable. Recent controlled truck-mounted experiments by Newton 
(1977) and Wang ^ (1980) agree closely with the previous 

findings. 

Schmugge (1980a) developed a relationship between soil texture 
and field capacity 

FC = 0.30 - 0.0023 x SAND + 0.005 x CLAY (22) 


where: 

FC = field capacity (for use with volumetric soil moisture) 
SAND = percent sand 
CLAY = percent clay. 
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He determined that passive microwave response should be related to 
percent field capacity instead of soil moisture by volume. Percent 
field capacity may, in some applications, be more important than 
absolute water content. 

The sensitivity of radiometers to soil moisture has also been 
demonstrated from space. McFarland (1976) showed a definite 
relationship between the Skylab 21 cm brightness temperatures and an 
antecedent precipitation index (API). The soil moisture index (API) 
was calculated frocn rain gage data within the 110 km footprint. 
Eagleman and Lin (1976) also related Skylab 21 cm data to soil 
moisture variations as determined by water budget models. The 1.55 cm 
Electrically Scanning Microwave Radiometer (ESMR) response was 
correlated to moisture indices (McFarland and Blanchard, 1977; Theis 
1979; Schmugge ^ a_[. , 1977). This shorter wavelength radiometer 
demonstrated the potential to detect soil moisture during periods of 
the year when the agricultural fields were essentially bare. 

Visible/Infrared Biomass Models 

Reflected solar radiation is the primary source of energy for 
visible/infrared sensors. Different materials possess different spec- 
tral reflective properties unique to that material. Due to Fresnel 
reflectance at air/water interfaces within leaves, radiation at the 
near- and middle-infrared frequency is strongly reflected (Gates, 
1980). This reflectance increases with increasing biomass. This 
phenomenon has been the basis for classification of vegetation and 
estimating biomass. Two biomass estimation models are presented here. 




Richardson and Wiegand (1977) developed the perpendicular vegeta- 
tion index (PVI). This graphical approach presented in Fig. 4 plots 
two Landsat bands against each other; MSS5(0. 6-0.7 i/n) and 
MSS7 (0.8-1. 1 pm). As soil reflectance properties change, the two 
bands change together to form the soil line. PVI was calculated using 
the Pythagorean theorem; 

PVI = \/(Rgg5 - Rp5)2 + (Rgg7 - Rp7)2 (23) 

where Rp is the reflectance values of an individual band and Rgg is 
the corresponding soil reflectance for the Landsat bands. PVI is the 
perpendicular distance from the soil background line to the plotted 
radiance values. 

The Transformed Vegetation Index (TVI) has been used as a range- 
land biomass estimator (Rouse et £l_. , 1973); 


TVI 



MSS7 - MSS5) 
MSS7 + MS$5!' 


+ 0.5 


(24) 


The ratioiiq under the radical is a standard practice to remove the 
scene to scene variations of sun angle and atmospheric optical density 
difference. The 0.5 was added to prevent the term under the square 
root from becoming negative. PVI was used in this study because it 
was developed for agricultural crop lands and may be used as both a 
crop classifier and biomass estimator. 


Sensor Co«bi nations 

Very few investigations have combined responses from various 
sensors to estimate soil moisture. Newton (1977) used the difference 
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FIG. 4 Diagram illustrating the principle of the perpendicular 

vegetation index (PVI) model. A perpendicular from candidate 
plant coordinates (Rp5, Rp7) intersects the soil background 
line at coordinates (RgS, Rg7). A PVI*0 indicates soil, and 
a PVI>0 indicates vegetation. From Richardson and Wiegand 
(1977). 
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in the horizontally and vertically polarized L-band radiometer at 35 
degree incident angle as a parameter to infer roughness. This 
parameter was used in combination with the L-band emissivity at 20 
degree incident angle as an estimate of soil moisture over various 
roughnesses. 

Jackson ^ (1981b) developed a model to predict emissivity of 

a scene from surface temperature, soil moisture, and vegetative bio- 
mass inputs. Encouraging results were obtained for wet biomass values 
below 1300 g/m^. Beyond this, the sensitivity is probably too low for 
the model to be of use. These results support the need for biomass 
estimations (possibly obtained from visible/infrared sensors) as an 
indicator of vegetation effect on bare soil emission. Jackson could 
not investigate this because visible/infra-ed data were not taken 
simultaneously with the microwave data. The data collected for this 
study were used to develop sensor combination techniques to estimate 
soil moisture by remote means. 

DATA COLLECTION AND PROCESSING 

Multisensor data over a wide range of frequencies are available 
only from NASA's C-130 aircraft. The C-130 collected all sensor data 
used in this study. Aircraft and soil moisture data were collected 
near Dalhart, Texas in August 1980 and nea." Guymon, Oklahoma in 
August, 1978. These areas were chosen because each had relatively 
homogeneous soils, reasonably low relief, a variety o' crops, and the 
availability of irrigation. The Texas and Oklahoma Panhandle areas 
were initially surveyed in one mile square sections with boundaries of 
the sections oriented North-South and East-West. Aircraft flight 
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lines were also orientated in these directions to make it possible for 
the pilot to fly one half mile over the fields that were sampled and 
thus allow adequate sampling by the sensors over each field. 
Agricultural practices usually dictate one crop per quarter section 
which is also the area irrigated by one center pivot sprinkler 
system. For this study a field was defined as one half of a quarter 
section with the long axis of the field parallel to the aircraft 
flight line. The half quarter section provided an area large enough 
for sufficient coverage with aircraft line sensors and small enough 
for an adequate ground sampling network. 

Experinent Site Descriptions 

Guymon 

SainplPd fields were located under four flight lines approximately 
2U km southwest of Guymon, Oklahoma (Fig. 5). The soil type was 

generally a silty clay (35% clay, 35% silt, and 30% sand) with many 
areas having caliche (CaCOj) near the surface. The crop types were 
alfalfa, and milo, with some fallow fields. The milo fields were 
chosen with some fields having rows perpendicular and some parallel to 
the flight lines. The field conditions are summarized in Table 1. 

Dalhart 

The sampled fields were chosen under two flight lines approxi- 
mately 20 km northwest of Dalhart, Texas (Fig. 6). The soil type of 
the surface 30 cm was a sandy loam (75% sand, 10% silt, and 15% 
clay). The 22 fields sampled were almost equally divided into 
irrigated and non-irrigated fields. The irrigated fields included 
millet and mature corn. Pasture, wheat stubble, and disked wheat 
stubble comprised the non-irrigated fields. 
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field maps, a) Area map of Guymon showing relative 
ns of each field map. 
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FIG. 5b Legend for Guymon, Oklahoma field maps. 




FIG. 5c Locations of sample 






ORiGlKAL PAGE 13 
OF PCH-'I? QUALITY 



sample fields at Guymon, 




ORIGINAL PAGE IS 
OF POOR QUALITY 


I 



33 


5e Locations of sample fields at Ruymon, 
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OF POOR QL'AUTi’ble 1. Surface Conditions of Guymon Fields 


Condition 


Fielcis Numbers 


Summer Fallow (bare) 

Sprinkler Irrigated August 1 
Sprinkler Irrigated August 14 
Listed August 17 

Circular Irrigated Alfalfa 


4, 13, 22, 27 


Circular Irrigated 
23 cm tall, drilled 
90 cm tall perper dicu'iar rows 

Furrow Irrigated 
90 cm tall parallel rows 
90 cm tall perpendicular rows 


7, 15 
20, 25 


8, IX, la, 2a 
18, 24 


Table 2. Surface Conditions of Dalhart Fields 


Condition 


Field Numbers 


Wheat stubble 
Disked stubble 

Mulched stubble (with weeds) 

Half disked, half burned stubble 


17,18 

19,20,21,22 

13,14 

15,16 

1,2,7,8,9,10,11,12 


Pasture 


Millet 
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DALHART, TEXAS 1980 

LEGEND FOR FIELD MAPS 1,2 & 3 


CROP 


Bare: wheat stubble 

disked wheat stubble 
mulched wheat stubble 

Millet 
Milo 

— Flight lirre markers 
A Corner reflectors 

# Rain gauges 

* Vegetation sample sites 
Row direction was east-west for all sample fields with row 


WF 

I I I I 





Corn 

Alfalfa 

Pasture 

Grazed 


crops. 




7 MILES 


by tK« Ttaot A&M UnivtrtiYy R*mot« S«niin 9 Bote compiUd from USCS to^ 09 repbic mops, 

R.S.C. toom fioM r^ofot on^ NASA cenfroctod •onel pKotogropby colUctod Awpwftt 14.18, 1980. 


FIG. 6b Legend for Dalhart, Texas field maps. 
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FIG. 6c Locations of sample fields at Palhart 
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FIG. 6d Locations of sample fields at Dalhart. 
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The summer of I960 was extremely dry In the Dalhart area. 

However, It was fortunate that this area received a one Inch rain 

after the first sampling date. This provided a range of soil moisture 
conditions. Those fields with vegetation were reasonably uniform In 
crop cover. All corn fields were In mature stages. Only fields 1 and 

2 showed areas of stress especially where some of the sprinkler 

nozzles were not functioning properly. Field conditions are 
summarized In Table 2. 


Soil Hoisture Data 

Soil moisture samples were collected at eight points per 32 
hectare field (Fig. 7) except In the Guymon fields 6 and 14 which were 
Intensively sampled (37 points). Oalhart sample depths were 0-2 cm, 
2-5 cm, 5-15 cm, 0-15 cm, 15-30 cm, 30-45 cm (Fig. 8). The 5-15 cm 
interval was divided Into two samples at Guymon; 5-9 cm and 9-15 cm. 
Manpowe- limitations did not allow each field to be sampled during the 

To correlate the ground data to aircraft data, the soil moisture 
values in fields not sampled on flight days had to be extrapolated. 
To do this, graphs were generated plotting soil moisture versus time. 
By incorporating precipitation events and irrigation schedules in the 
graphs, it was possible to produce dry-down curves for each field and 
point. This method also enabled the estimation of soil moisture for 
points where values were missing due to pesticide spraying, rain, 
sample spills, etc. Where the curve for each depth intersected a 
flight time, the value was recorded as the new gravimetric soil 
moisture. These values plus the bulk density results were used to 
calculate the volumetric soil moistures by print. Texture analyses 
were performed for each field In Guymon and for selected fields In 


Oalhart. 






805 METERS 


These points weie moved outside the pivot houndary lor non>citculaf fields. 


FIG. 7 Sampling pattern for fields at liuymon and Dalhart. 


41 


original page b 





Ground tracks of the line sensors were used to determine which 


sample points fell within the area covered by the sensors. These 
points were arithmetically averaged by field and termed soil moisture 
under the line sensors or just soil moisture. 

The soil sample processing lab was set up with conventional 
microwave ovens. As the samples came in from the field, they were 
weighed with seals and lids on the cups. This weight (measured to a 
tenth of a gram) was recorded as the wet weight. Lids and seals were 
then discarded and soil samples dried in the ovens. At regular time 
intervals two test samples per oven were removed and weighed and 
returned to the ovens. The weights were plotted on a graph to track 
the dry-down. When the weight loss of the test samples leveled out to 
less than a 0.2 gram difference frun the previous weighing, they were 
considered dry. All samples were then removed from the oven and their 
dry weights recorded. As the preprocessing was completed for a set, 
the samples were stored for possible future reference. 

A programmable calculator was used in the lab to compute gravi- 
metric soil moistures as the final weights were completed. The pro- 
gram automatically subtracted an average lid, seal and cup weight from 
the gross values. Daily soil moisture charts were kept on each field, 
point and depth. This system facilitated quality control of the 
data. If a value was questionable, the sample was located, reweighed 
and its soil moisture recalculated immediately. 

Aircraft Data 

Aircraft sensor data were similar for both experiments. The 
Guymon data were collected on August 2, 5, 8, 11, 14, and 17, 1978. 
The Dalhart data were collected on August 14, 16 (two flights), and 18 
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(late afternoon), 1980. The sensor package consisted of 1) seven 
scatterometer frequencies and polarizations, 2) three passive 
microwave radiometer frequencies and polarizations, 3) Barnes PRT-5 
thermal data, 4) aerial photography, and 5) visible/infrared scanner 
data. Soil moisture flights were made at an altitude of 500 m. 

Scatterometer frequencies and polarizations included: 


1) 

13.3 

GHz 

VV 

2) 

4.75 

GHz 

HH 

3) 

4.75 

GHz 

HV 

4) 

1.6 

GHz 

HH 

5) 

1.6 

GHz 

HV 

6) 

0.4 

GHz 

HH 

7) 

0.4 

GHz 

HV 


where VV = vertically polarized transmitted and received 
HH = horizontally polarized transmitted and received 
HV = horizontal !> polarized transmitted and vertically received. 
Eight look angles from nadir were processed: 5, 10, 15, 20, 25, 35, 
40, and 45 degrees. Near-nadir (3 degree) passive microwave data were 
collected at 1) L-band (1.6 GHz) horizontal polarization, 2) C-band 
(4.75 GHz) vertical, and 3) C-band horizontal polarization. In 
Dalhart, the eight channels of the NSOOl (simulated thematic mapper) 
bands included: 

1) Channel 1: 0.45-0.52 ym 

2) Channel 2: 0.52-0.60 gm 

3) Channel 3: 0.63-0-69 ym 

4) Channel 4: 0.76-0.90 

5) Channel 5: 1.00-1.30 ym 

6) Channel 6: 1.55-1.75 ym 


44 


7 ) 

8 ) 

In Guymon, five channels 
were processed: 


Channel 7: 2.08-2.35 
Channel 8: 10.4-12.5 pm 
from the modular multispectral 


scanner (M^S) 


1) Channel 4 0.548-0.583 pm 

2) Channel 7 0.662-0.701 pm 

3) Channel 8 0.703-0.747 pm 

4) Channel 9 0.770-0.863 pm 

5) Channel 11 8.000-12.00 pm 

The M^S channels 7, 9, and 11 correspond closely to the NSOOl channels 
3, 4, and 8, respectively. 


Data Processing 

Processing procedures for the individual sensors and ground data 
are given in Appendix A. The procedures included calibration, quality 
control, data omission, and field average determination. Data omis- 
sion usually occurred with excessive roll and drift of the aircraft. 
Once field averages had been calculated for each sensor and their cor- 
responding soil moisture, the ground and aircraft data sets were 
merged. Further analysis was done on the Field averaged data set. 

ANALYSIS AND DISCUSSION 
Approach 

Case studies were first performed to investigate the microwave 
sensor responses to soil moisture and roughness. Previous studies 
have related microwave sensor responses to surface soil moisture at 
near incident angles. Near incident angle scatterometer (active) and 
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radiometer (passive) microwave data for all available frequencies and 
polarizations were plotted and regressed against soil moisture (0-2 
cm) over essentially bare soils. Regression line slopes and 
correlation coefficients were used as criteria for evaluating the 
sensitivity and relative value of each sensor configuration. From 
this analysis a single microwave sensor, the L-band radiometer, was 
chosen as the sensor that exhibited the highest potential to quantify 
soil moisture. 

The analysis was then expanded to vegetated fields for this 
sensor. Linear relations were determined for the L-band response to 
soil moisture for each crop type. These individual relations were 
used to obtain predicted values of soil moisture. The correlation 

between measured and predicted soil moisture over all fields 
represented the optinum capability of the L-band passive microwave 
sensor to quantify soil moisture when crops can be accurately 
clafsified. 

In general, crop type is related to biomass; however, there is a 
range of biomass within each crop type. A visible/infrared biomass 
estimator and crop classifier (PVl) was directly combined with the 
L-band passive microwave data to quantify soil moisture independent of 
crop type. The direct combination technique was comfiared to the 
classif ication techniques and microwave alone. 

Case Studies 

Two case studies were performed on the sensor data from Dalhart 
in order to Investigate each microwave sensor's response to particular 
changes In 1) soil moisture and 2) roughness. Fig. 9 Illustrates a 







field being actively irrigated at the time of the data collection. 
Soil moisture measurements were not available but the field is notice- 
ably wetter (darker) on the left side where it had been irrigated 
recently. The field had unifonn roughness on a scale similar to near- 
by sample fields. Line plots of each sensor are presented at the same 
horizontal scale as the photograph in Figs. 10-12. Passive microwave 
sensor (Fig. 10) exhibited approximately a 25 K sensitivity to the 
difference in soil moisture while the surface thermal temperature 
range was approximately 3.5 K. This small thermal response is due to 
the time of day the data were taken. The flight time was mid-morning 
but this change in thermal temperature indicates that there was a 
change in moisture between the two portions of the flight line. The 
scatterometers at 10 degree incident angle (Fig. 11) exhibited con- 
siderable signal noise and a maximum sensitivity of only approximately 
5 dB after the noise is smoothed. It is interesting to observe that 
all the scatterometers exhibit a stronger dependence on soil moisture 
at 40 degree incident angle (Fig. 12). This is contrary to results 
from most previous studies. 

The second case study involved field 15 at Dalhart. A photo is 
not presented but the left half portion of this field had been 
recently disked. The right half was burned wheat stubble (relatively 
smooth). The field had unifonn soil moisture (15 percent by volume). 
A definite roughness influence i". apparent in the passive microwave 
responses (Fig. 13) that does not readily appear in the 10 degree 
scatterometer data (Fig. 14). The rojghness effect is again apparent 
in all the 40 degree scatterometer data (Figure 15) with the exception 
of returns from the 1.6 GHz HV, 0.4 GHz HH, and 0.4 GHz HV scattero- 
meters. 
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FIG. 11a Line plots (sensor response vs. time) of field in Fig. 9 for 
13.3 GHz VV, 4.75 GHz HH, and 4.75 GHz HV scatterometers at 
10 degree incident angle. 
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FIG. 12a Line plots (sensor response vs. time) of field In Fig. 9 for 
13.3 GHz VV, 4.75 GHz HH, and 4.75 GHz HV scatterometers at 
40 degree incident angle. 
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FIG. 12b Line plots (sensor response vs. time) of field in Fig. 9 for ^ 

1.6 6H2 HH, 1.6 GHz HV. 0.4 GHz HH, 0.4 GHz HV scatterometers I 

at 40 degree incident angle. 
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FIG. n Line plots (sensor response vs. time) of case study field 
(rough and smooth) for L- and C-band radiometers and PRT-5. 
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nu. 14a Line plots (sensor response vs. time) of case study field 

(rough and smooth) for 13.3 GHz VV, 4.75 GHz HH, and 4.75 GHz 
HV scatterometer at 10 degree Incident angle. 
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FIG. 14b Line plots (sensor response vs. time) of case study field 
(rough and smooth) for 1.6 GHz HH, 1.6 GHz HV, 0.4 GHz HH. 
and 0.4 GHz HV, scatterometer at 10 degree incident angle. 
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FIG. 15b Line plots (sensor response vs. time) of case study field 
(rough and smooth) for 1.6 GHz HH, 1.6 GHz HV. 0.4 GHz HH, 
and 0.4 UHz HV scatterometer at 40 degree Incident angle. 





Microwave Sensor Responses to Soil Moisture Over Bare Fields 

This analysis is presented in two sections (scatterometers and 
radiometers) to avoid confusion between the active and passive sys- 
tems. For the same set of fields, plots of sensor responses versus 
volumetric soil moisture (0-2 cm) are presented for each sensor. The 
fields selected for these illustrations were all bare with a range of 
roughness effect due to tillage and surface weathering. There was no 
quantitative measure of roughness, but it was observed to range from 
smooth to moderately rough as defined by Newton (1977). 

Scatterometers 

As mentioned before, previous investigations indicated that the 
scatterometer's highest potential for detection of soil moisture is 
around the 10 degree incident angle. Plots of each scatterometer's 
response at the 10 degree incident angle versus soil moisture in the 
surface 2 centimeters are presented (Fig. 16-19) to illustrate 
differences in the data sets. The scatterometer returns at low soil 
moisture shifted approximately 0 dB, 5 dB, 8 dB and 7 d3 for the 13.3, 
4.75, 1.6 and 0.4 GHz systems, respectively. These shifts illustrate 
that the comparison of the two data sets must be approached with 
caution. Further analysis was conducted for each scatterometer within 
each data set to minimize the influence of lack of calibration. 

Each scatterometer's response versus volumetric soil moisture 
(0-2 cm) was analyzed for 10, 15, 20, and 40 degree incident angles. 
The 1.6 GHz HH scatterometer (Dal hart) is illustrated in Fig. 20 as a 
representative sample of the scatterometer responses. The remainder 
of the illustrations are in Appendix B. The absence of a direct 
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FIG. 19 Scatterplot of 0.4 GHz o° vs. volumetric soil moisture (0-2 cm) for Guymon 
and Dalhart bare fields. 
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FIG. 20 Oalhart 1.6 GHz HH o° vs. volumetric soil moisture (0-2 cm) 
at 10, 15, 20, and 40 degree Incident angles. 
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and soil moisture for all bare fields and each scatterometer (10, 15, 
20 and 40 degree incident angles) are compiled in Table 3. 


Table 3. Scatterometer Values for o° versus 
Volumetric Soil Moisture 



R^ Values 


Guymon 

Dal hart 

Scatterometer 

10° 

15° 

O 

O 

OsJ 

0 

O 

10° 

15° 

ro 

o 

o 

C 

O 

13.3 VV 

0.43 

0.47 

0.65 

0.10 

0.20 

0.29 

0.21 

0.22 

4.75 HH 

0.32 

0.35 

0.46 

0.60 

0.13 

0.14 

0.12 

0.03 

4.75 HV 

0.57 

0.36 

0.49 

0.51 

0.45 

0.49 

0.45 

0.36 

1.6 HH 

0.48 

0.58 

0.47 

0.35 

0.27 

0.46 

0.43 

0.28 

1.6 HV 

0.43 

0.59 

0.57 

0.62 

0.02 

0.28 

0.36 

0.17 

0.4 HH 

0.02 

0.11 

0.19 

0.22 

0.08 

0.38 

0.19 

0.07 

0.4 HV 

0.02 

0.00 

0.00 

0.00 

0.29 

0.49 

0.22 

0.17 


* The largest value (0.65) obtained for either data set was 

Guymon's 13.3 GHz VV (20 degree); however, the largest average 
value (0.52) was obtained for the 1.6 GHz HH (15 degree). These 
f values will be used in the next section to compare the capability of 

the scatteroineters to detect soil moisture to that of passive 
microwave radiometers over the same set of fields. 


Radiometers 

Radiometers provide intensity measurements in terms of brightness 
temperature (K). By assuming the temperature of the emitting layer to 
be that sensed by the PRT-5, emissivity can he calculated using 
equation (18) as 

£ = Tbt/TprT-5 ( 25 ) 


1 
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For this study the emissivity was used as the measure of the passive 
microwave response. The data sets were Initially analyzed together 

(Figs. 22-24) to determine if there were calibration problems. No 
significant calibration problems were evident except In the L-band 
horizontal emissivity (Fig. 24) where the dry emissivity values were 
significantly higher at Dalhart; however, the linear value for this 
combined data set was 0.71. The radiometer data were further analyzed 
by comparisons to volumetric soil moisture and percent field capacity 
for each data set (Figs. 25-27). Both the C-band horizontal and 
vertical data showed considerable scatter for Dalhart (Figs. 25 and 
26). The scatter is related to the differences in the roughness of 
the field groups. Fields D17&18 appeared "smoother" to the 
radiometers tnan fields D15&16. Recall that these two fields also 
exhibited corresponding differences In relative roughness In the 40 
degree Incident angle 1.6 GHz HH scatterometer (Fig. 20). The same 
roughness effect is apparent In the L-band radiometer data but to a 
lesser degree (Fig. 27). The "sensed" relative roughness did not 
Influence the longer L-band wavelengths as much as it did the C-band 
wavelengths. 

Percent field capacity (PFC) is defined as 

PFC = VSM/FC X IJO (26) 


where 

VSM = measured volumetric soil moisture 

FC = Field capacity obtained from texture analysis and equation 

( 22 ) 
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FIG. 22 Scatterplot and regression line for C-band (vertical 

polarization) emlsslvlty vs. volumetric soil moisture for 
Dal hart and Guymon bare fields analyzed together. 
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FIG. 23 Scatterplot and regression line for C-band (horizontal 

polarization) emissivity vs. volumetric soil moisture for 
Dal hart and Guymon bare fields analyzed together. 
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FIG. 24 


Scatterplot and regression line for L-band (horizontal 
polarization) emlsslvlty vs. volumetric soil tftolsture for 
Dalhart and Qiiymon bare fields analyzed together. 
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FIG, 25 C-band (vertical polarization) response to volumetric soil 
moisture and percent field capacity for Guymon and Dalhart 
bare fields analyzed separately. 
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26 C-band (horizontal polarization) response to volvnetric soil 
moisture and percent field capacity for Guymon and Dalhart 
bare fields analyzed separately. 
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FIG. 27 L-band (horizontal polarization) response to volumetric soil 
moisture and percent field capacity for Guymon and Dalhart 
bare fields analyzed separately. 



The difference in soil texture betvyeen the Guymon and Dalhart test 
sites is better identified in the L-band radiometer measurements since 
the effect of roughness, which could be confused with texture effects 
is reduced at this frequency. Guymon had a much higher clay content 
than Dalhart. The slopes of the regression lines are different at 
Guymon and Dalhart when e is related to volumetric soil moisture 
(-0.0074 and -0.0011, respectively). In contrast, they are almost 
identical when e is related to percent field capacity. This supports 
the hypothesis that passive microwave responses should be analyzed as 
a functio". of percent field capacity. After consideration of the 
field capacity plots it is evident that the calibration problem with 
L-band was magnified in the percent field capacity plots because it 
was not offset by the texture effect. The difference in the intercept 
between data sets is 0.043. This correction was added to the Guymon 
data set for all further analyses. The correction was added to Guymon 
because the L-band antenna used in collecting that data set was 
deteriorating and replaced between the times of the Guymon and Dalhart 
experiments. The deterioration tended to reduce the passive microwave 
response at Guymon (Blanchard, 1981). 

The one inch rainfall which occurred at Dalhart provided an 
opportunity to relate the roughness "sensed" by L-band radiometers 
with that sensed by L-band scatterometers. This was not possible in 
Guymon because the irrigated bare fields (G6&14) provided only one 
field with a sufficient moisture range to establish the relative 
roughness effects. The o° intercept values at 0 percent moisture, 
0 °', were obtained for each of the Dalhart bare field groups from the 
1.6 GHz HH (40 degree) scatterometer (Table 4). 


74 


Ch«Gl*\AL 

OF POOR QUALITY 


Table 4. oq' Values for Each Field Group in Dalhart 


Field 


013&14 

ni5S16 

D17S18 

019420 

020421 


-30.8 

-28.7 

-36.2 

-32.6 

-31.5 


This "roughness" was considered to influence both the intercept and 
the slope of the c vs percent field capacity relationship. The linear 
equation: 


was modified by 


e = A + B * PFC 


A = A1 i A2 * Oq' 


(27) 


(28) 


B = B1 ♦ B2 * Oq‘ 


(29) 


c * A1 ♦ (81 8 PRF) + (B2 * PFC * ) * (A2 * o^' ) (30) 


The resulting equation when applied to bare fields in the Dalhart data 
set is: 

e « 1.212 ♦ (0.00165*PFC) ♦ (0.l)00145*PFC*o°' ) ♦ (.00681 x o°') (31) 

with value = 0.93 and F value = 158.5. The significant increase in 
R^ values indicates that scatter in the L-band radiometer is indeed 
due to difference in relative roughness as sensed by the 1.6 HH (40 
degree) scatterometer. From the case study and comparison of 
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values, the L-band radiometer shows the strongest relationship to soil 
moisture over the same set of bare fields than any other microwave 
sensor. Therefore, L-band radiometer data were used in all further 
analyses. 

L-band Emissivity Response to Soil Moisture over all Crop Types 

No significant difference in e response to moisture could be dis- 
cerned when certain "near bare" fields were included in the classifi- 
cation of bare fields. These included Guymon's 23 cm tall milo fields 
(7 and 15) and Dalhart's pasture (5 and 6) and millet (3 and 4). This 
indicates that for L-band, attenuation of the soil emission is insig- 
nificant for these fields. For this experiment the crop types were 
classified and generally can be rank-ordered from lower to higher 
biomass: 

1) bare soi 1 

2) alfalfa 

3) milo 

4) corn 

Four fields of alfalfa were sa'itpled at Guymon. Fields 4 and 13 
and fields 22 and 27 were pairs of adjoining fields. During the 18 

day experiment the alfalfa was irrigated and matured. Fields 22 and 
27 were harvested before the last flight day. These two fields were 
on a steep hillside; consequently, the effective incident angle was 
increased from 3 degrees to greater thon 10 degrees. An increase in 
effective incident angle in the horizontal polarization causes a 
decrease in emissivity (Newton, 1977). The decrease was observed in 
Fig. 28a where the values of c for fields G22427 were generally 
shifted downward from the values for fields G4A13. A realistic 
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correction of 0.02 (Newton, 1981) was added to the emissivity for 
fields G22&27 with a resulting increase in the significance of the 
linear relationship between e and PFC (Fig. 28b) (R^ values increased 
from 0.49 to 0.67), All milo fields considered in this study were at 
the Guymon site while all corn fields were at the Dal hart site. The 
linear analysis of the milo fields (Fig. 29) showed a strong relation- 
ship between c and PFC (R^ value = 0.69). The same analysis for corn 
fields (Fig. 30) yielded a much weaker relationship (R^ value = 0.08) 
due to: 

1) Corn masked or attenuated the emissivity from the soil, thus 
destroying the response to soil moisture differences. 

2) The difference of biomass depending on the field group varies 
markedly within the corn classification. 

The influence of different biomass within the corn fields is 
noted when PVI values are compared to regression lines for individual 
field groups. Fields 07&S and 011&12 had PVI values greater than 
4.5. These fields were combined into field D>4.5. Regression lines 
for the three subclasses of corn are presented in Fig. 31. The 
relatively flat slope for the combined data for high biomass corn 
fields (D>4.5) indicates thet as PVI values approach 4.5 the measured 
emissivity is virtually insensitive to the soil moisture under the 
corn. To illustrate the need to compensate for vegetation, emissivity 
was plotted versus percent field capacity for all fields combined with 
the resulting R^ values of O.Oh (Fig. 32). 

Most previous studies have modeled emissivity as the dependent 
variable while using soil moisture and biomass as independent 
variables (Jackson 1981b; Choudhury, 1981). Choudhury used 

the model 
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FIG- 29 ScaUerplot and regression line for L-band emissivlty vs. 
percent field capacity for Milo fields. 
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FIG. 30 Scatterplot and regression line for L-band emlsslvity vs. 
percent field capacity for Corn fields. 
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FIG. 31 Scatterplot and regression lines with corresponding PVI 

values for L-band emissivity vs. percent field capacity for 
individual fields of corn. 
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FIG. 32 Scatterplot of L-band emissivlty vs. percent field capacity 
for an crop types. 
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G^e'^ ♦ (1 - Gj) Cy(l - e'^) e' 


♦ Cy(l - e"^) 


( 32 ) 


where t = vegetation Index 

Gm ' measured emissivlty 

Gs » emissivlty of the soil (function of soil moisture) 

Gy “ emissivlty of vegetation 

The first tenn on the right Is the contribution from the soil surface 
as attenuated by the vegetation. The second term Is the contribution 
from the vegetation as reflected by the soil surface (1-es) 
again attenuated by the vegetation. The third term Is the direct 

contribution from the vegetation. By rearranging terms: 


^m 


- Ey 


c, e- . G, e-2^) 


( 33 ) 


and solving for €5. the resulting equation is: 


e 


s 




c 


e 


-2t 

V 



( 34 ) 


There Is always an error In the measurement of Gm. cm was 
replaced with g^ + G^r (ter ’S the error in measurement) 

yielding 


G 


s 





-2^ 






( 35 ) 


In a similar fashion, the error term can be applied to the emissivlty 
of the soil Gs resulting in the following equation 
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e, (c^ . Cg^)(e- 




(36) 


The soil moisture information is contained in cs* Choudhury 

reported values of cy to be 0,97. Using the value eer * 0.01, the 

sensitivity of cs and ep to different values of t was investigated 
in Fig. 33. It was noted that as t increases from 0 to 1.5 the error 

in es increases from 0.01 to 0.18 while the error of cm decreases 

from 0.01 to 0.0006, This indicates that these equations can predict 
the cmissivity measured at the radiometer well, but cannot estimate 
the e$ or emissivity of the soil accurately when the vegetation 
becomes dense. 

Coefficients of the linear regression equations for the calcula- 
tion of e from percenc field capacity for each crop type are tabulated 
in Table 5. 

Table 5. Coefficients of linear prediction equations for 
estimation of emissivity as a function o^ percent 
field capacity. 


Crop Type 

Intercept 

Slope 

Bare 

0.98U1 

-.00253 

Alfalfa 

0.9638 

-.00125 

Milo 

0.9069 

-.00136 

Corn 

0. 9657 

-.000121 


These coefficients were used to predict the value of emissivity. The 
calculated values of emissivity resulting from these equations are 
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FIG. 33 Illustration of the sensitivity of Choudhury's model to 

errors in measurement vs. a vegetation index (t). The curve 
denoted by cj is the error in the predicted emissivity of 
the soil that would be incurred with a 0.01 error in the 
measured emissivity. The curve denoted by tn, is the error 
in the predicted measured emissivity with a 0,01 error in 
cs* 





compared to the measured emissivity in Fig. 34. The value of 
measured versus predicted emissivity using this approach is 0.74. 

Soil Holsture Determination with Multiple Sensors 

The remainder of this study considers the soil moisture parameter 
as the dependent variable. Regression lines for percent field capac- 
ity versus emissivity for the four major crop types along with mean 
vegetation index values are shown in Fig. 35. This figure introduces 
two approaches to estimate soil moisture by using visible/infrared 
data with passive microwave sensors. The first approach uses 
visible/infrared sensors to classify the crop and crop class to esti- 
mati moisture. The second approach bypasses cl assi 1 i cation by using 
the PVI values (obtained from visible/infrared sensor) and microwave 
response in a direct relationship. 

Classification Technique 

The classification technique approach assumes that all fields 
within a particular crop type have equal biomass and thus the same 
vegetation effect on the passive microwave radiometer. This assump- 
tion was tested by investigating the effect of different biomass (as 
indicated by PVI) within crop types for corn and milo. Regression 
lines for each field group with the corresponding average PVI values 
are presented in Fig. 36. Significantly different regression lines 
within both crop types were observed. The sensitivity of the percent 
field capacity to changes in emissivity was directly related to PVI. 
That is, as PVI increased, the sensitivity of PFC to c increased. 
This increased sensitivity dictates that errors in the measurement of 
e have greater effects on the accuracy of estimates of the percent 
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FIG. 36 Regression lines relating percent field capacity to emissivity and the PVI values 
for each field group in a) milo and b) corn. 


field capacity at high PVI values. This same relation is the control- 
ling element in the sensitivity analysis of the Choudhury equation 
that was discussed in the previous section. 

Coefficients for the linear regression equations used for esti- 
mating percent field capacity as a function of measured emissivity are 
presented in Table 6. The regression lines defined by these values 
were illustrated in Fig. 35. 

Table 6: Coefficients of linear prediction equations for 

estimation of percent field capacity as a 
function of emissivity 



Crop Type 

Intercept 

Slope 


Bare 

291.86 

-291.97 


Alfalfa 

493.61 

-493.65 


Mi lo 

512.70 

-510.19 


Corn 

707.31 

-656.58 


Measured percent field capacity was plotted against the predicted 
values obtained from these equations (Fig. 37a). This illustration 
points out the failure of using one set of linear coefficients per 
crop type (especially in corn) to predict soil moisture for an entire 
set of fields inside that crop type. The most significant value 
which can be expected using classification technique on this data set 
with corn included is 0.67. When the corn fields were not considered 
the value increased to 0.76 (Figure 37b). Thus, if the severe 
attenuation from corn can be avoided we may be able to estimate per- 
cent field capacity from measured emissivity reasonably well. 
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Direct Combination 


The estimation of percent field capacity with a combination of 
sensors required the development of an equation that defines percent 
field capacity as a function of both microwave emissivity and vegeta- 
tion index data (PVI). This technique should be more flexible than 
the classification technique described earlier since it can accommo- 
date differences in biomass within each crop class. 

The regression lines for the different crop types and the corres- 
ponding values of PVI are presented in Fig. 38 with the corn sub- 
divided into three fielo groups. The slopes of these regression lines 
were then plotted against the average PVI values in Fig. 39. The 

form of the relation appears to be a straight line function when 

fields D9&10 and D>4.5 are not considered. The linear form of the 
relationship between percent field capacity and emissivity is: 

PFC = A - B * c (37) 

If biomass (PVI) is considered to linearly affect both A and 8 then: 

A = A1 + A2 * PVI (38) 

B = B1 + B2 * PVI (39) 

Substituting into equation (37) and expanding the form, the equation 
becomes: 

PFC = A1 + (A2 * PVI) - (B1 * c) - (B2 * e * PVI) (40) 

Applying this to all fields except those with average PVI value 

greater than 4.3 (D9&10 and 0>4.5) yields the equation: 

PFC = 279.53 + (51.20 * PVI) - (281.22 * e) - (48.41 * c * PVI) (41) 
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FIG. 38 Regression lines relating percent field capacity to 

emissivlty for all fields (corn divided Into Individual 
field groups). 





It is significant to note the values of A1 ani B1 in equation 
(41) are comparable to the values of A end B obtained for bare fields 
(291.86 and 291.97, respectively) (Table 6). This indicates that the 
influence of vegetation was compensated for by the Increased values of 
PVI obtained over the non-bare fields. Comparisons of the capability 
of the L-band radiometer to estimate percent field capacity are 
presented in Fig. 40 when it is used: 

1) alone 

2) with a classification technique which excludes all corn 

3) in direct combination with PVI excluding fields with average 
values of PVI greater than 4.3. 

This illustrates that passive microwave alone cannot estimate soil 
moisture when several crop types are considered. The classification 
technique works reasonably well when the variation of biomass within a 
crop type are small. This technique also requires perfect classifica- 
tion. The direct combination technique shows the greatest potential 
to detect soil moisture because it: 

1) is able to discern difference in biomass within crop types 

2) will vary as the crop matures 

3) is not hindered by misclassification 

4) provides a means to determine a threshold value of biomass 
beyuiid >-hich the radiometer is unable to accurately estimate 
soil moisture. 
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CONCLUSIONS AND RECONMENDATIONS 


The principal objective of this research was twofold. The first 
objective was to select the most useful single sensor from the arr&y 
of sensors available for the estimation of soil moisture. The second 
objective was the development of a linear combination of senso'' data 
that would provide an Improvement In estimation of soil moisture 0”er 
estimates made with a single sensor. 

The L-band passive microwave radiometer was found to be the most 
suitable single system. Several subordinate conclusions led to the 
selection of this sensor. First the L-band system was found to be the 
most sensitive to soil moisture differences In terms of the range of 
return frotn the full range of soil moisture. Secondly, the error In 
measurement in relation to the range of sense' response was small In 
relation to the active microwave systems. Third, It became obvious 
that with the passive microwave systems It was possible to get a rela- 
tive calibration while with the active microwave systems this could 
not be done. Fourth, the L-band horizontal passive microwave res- 
ponded less to variations In roughness than the C-band passive micro- 
wave. In the preliminary analysis leading to these conclusions. It 
was shown that the L-band passive microwave measured emissivlty rela- 
tionship to soil moisture was more significant when moisture was 
expressed as percent of field capacity. 

All previous experiments analyzed the microwave return, whether 
it was the passive microwave antenna temperature or the active micro- 
wave scattering coefficient, as a function of the moisture measure- 
ments. A more useful but more difficult modeling effort is required 
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to predict moisture frotn remote sensing inputs. Due to the nature of 
the influence of surface roughness and vegetation on the transmission 
of microwave energy from the surface soil volume to the antenna, sig- 
nificant errors in estimated soil moisture will occur unless the sys- 
tems used can provide a reliable measure of surface roughness and 
vegetation attenuation. No practical satisfactory measure of rough- 
ness was found in these data; however, the perpendicular vegetation 
index was found to be useful as a measure of vegetation effects for 
agricultural crops other than dense corn. A linear equation was 
developed to predict percent field capacity as a function of L-band 
emissivity and the vegetation index. The prediction algorithm 

improved the estimation of moisture significantly over predictions 
from L-band emissivity alone. 

The results of this experiment infer that effective estimation of 
soil moisture will require as a minimum a combination of a dual polar- 
ized L-band radiometer {at some incident angle greater than 15 
degrees) and a visible/infrared imager with at least two bands that 
are suitable for developing the npppendicular vegetation index. Most 
importantly all of the sensors ust have a satisfactory means for 
calibration. An active microwave system could also be used for the 
roughness measurement if such a system is calibrated. Further experi- 
ments along these lines should be conducted with truck systems in 
order to maintain more control of the site and better calibration of 
instruments, while collecting and processing the data at lower cost. 
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DATA PROCESSING PROCEDURES 


The processing of the Guymon and Oalhart data sets was 
accomplished through joint effort at Texas ASM University. This 
appendix summarizes the work and can also be found In Blanchard 

(1980), McFarland (1981), and Rosenthal (1981). 

Scatteroneter Processing 

Scatterometer data were collected aboard the NASA C-130 In analog 
form on a 14-track tape. Copies of the tape were later sent to Texas 
ARM University/Remote Sensing Center for Individual processing, which 
consisted of three phases. The Initial processing converted the 
analog data to digital and copied the digital data onto 9-track 

magnetic tapes. The second phase processed the digital data using 

software which calculated sigma (o°) values for each look angle at 

given time intervals. Data were processed so that a cell size roughly 
had a length of 25 m for 13.3 GHz, 38 m for 4.75 GHz, 50 m for 1.6 
GHz, and 75 m for 0.4 GHz. The processing software was described by 
Claassen^^. (1979) and Clark and Newton (1979). 

After processing scatterometer data, field start and stop times 
were determined for each frequency and polarization from line plots of 
0 ° versus time, and aerial photographs. Times were adjusted by shift- 
ing the start/stop times at least 0.5 seconds to Insure full scattero- 
meter coverage within the field. The final start and stop times 
defined the field boundary and were used In determining field averages 
for each frequency, polarization, and look angle. Time frames during 
excessive aircraft roll and drift (roll greater than 3.5*; drift 
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greater than 9®) were noted and data from affected look angles were 
deleted from further analysis. 

No known technique or mechanism was available to calibrate all of 
the scatterometers. Consequently, any temporal variation In o* could 
Indicate either soil moisture, roughness, or vegetation changes, or 
unstable sensors. 

Due to excessive aircraft roll and drift, several look angles had 
to be eliminated at Dal hart and Guymon due to the uncertainty of the 
cell being within the field. At Dalhart, data from only one field had 
to be el1m1nated--f1eld 16 on 8/18/80. Also, data at 40® and 45® look 
angles off nadir from several other fields on 8/18/80 were eliminated 
due to excessive drift (Table Al). At Guymon, flying conditions were 
much worse; consequently, data from more fields needed to be deleted. 
A complete list of omitted look angles Is given In Table A2. Data 
from 8/11, 8/14, and 8/17/78 were most questionable. 

Passive Microwave Processing 

The raw analog data collected aboard the aircraft were converted 
to digital uncorrected brightness temperatures at NASA/Goddard Space 
Fight Center (GSFC). Corrected brightness temperatures were 
calculated from an equation developed at NASA/JSC (O'Neill, 1981): 

1 L r2(T )(L) 

'b ■ T - ■ r° — - ’l»-» - <'“> 

where t Is the transmittance of the radome, e Is the emissivlty of the 
radome, T(j Is the uncorrected brightness temperature based on raw 
digital counts, L Is antenna cable loss factor, T|^ Is an antenna 
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TABLE A!. QUESTIONABLE SCATTE«OK£TEK OATA FOR OALHART 
AS OF 5-6-81 


Date 


Field 

# 


Questionable Analysis 

8/14/80 

All 

data 

Is good 



8/16/80 

All 

data 

Is good 



8/18/80 

L12 

R2 

20,8,18 

45* 

(drift 9*) 


L12 

R2 

14 

40, 

45“ (drift 11“) 


Lll 

R3 

16 

All 

Angles 
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1 

I 


8/2/78 


8/8/78 


8/11/78 


8/14/78 


8/17/78 


LI 

R1 

2,4,6, 7.e.2x. lx 

L2 

R1 

10,13,14, 15,2a.2x,l 

LI 

R2 

2, 4, 6, 7, la, 2x, lx 

L2 

R2 

15,17,2a 

L2 

R1 

17, lx 

L2 

R2 

2A 

L4 

R1 

26 

LI 

R2 

2,6.7 

LI 

R1 

6,8,2x 

L3 

R1 

19,22, lx 

L2 

R1 

2x, 

L4 

R1 

24,25,27 

LI 

R2 

4,6,7,1A 

L3 

R2 

22 

L2 

R2 

10,17 



2A, 2X 

L4 

R2 

24,26,27 

LI 

R2 

4 

L3 

R2 

19 

L2 

R2 

13 



10 

LI 

R3 

all fields 

L3 

R3 

lx 

L2 

R3 

13,14 



15 

L3 

R1 

21,22 

LA 

R1 

2x, 24, 25,26, 27 

L3 

R2 

21,22 



lx, 19, 20 

L4 

R2 

24,25,2x 


R 


EROMETER DATA FOR GUYMON 


Questionable Analysis 


40* .45* (-8* drift. 2* roll) 
45* (-9* drift) 

45* (-9* drift) 

45* (-8* drift) 

all angles 
all angles 
all angles 
all angles 

all angles 
all angles 
all angles 
all angles 
all angles 

all angles 

45* (-4* drift. 4* roll) 
ail angles 
all angles 

all angles 

40*. 45" (-8* drift, 3* roll) 
45" (9" drift) 

40*. 45" (9" drift, 3" roll) 
40", 45" (11" drift) 
all angles 
all angles 
45" (9" drift) 

35*. 40" ,45* (-12* drift) 

35*. 40", 45" (-12" drift) 
all angles 

40" ,45" (-10" drift) 

45" (-9" drift) 





These fields were deleted from the MFMR plots due to excvssl.e roll; 
drift was not a factor. 

temperature factor, Tr Is the radome temperature factor, Is af* 
Internal parameter for each frequency, and U the self-emission 

of the receiver. The various constants used In the equation were 
determined from flights over homogeneous areas. Once brightness tem- 
peratures were calculated, line plots of Tg versus time were pro- 
duced and field start and stop times were determined from the .lots. 
The times defined field boundaries and were used In calculating field 
averages for each frequency and polarization. 

Since the passive microwave radiometer was oriented at a constant 
angle (3" from nadir), any excessive roll would Imply questionable 

MFMR data. Consequently, any time the airplane had roll greater than 

3.5®, the field average MFMR data we.*e deleted. Table A3 lists the 

deleted data. With the exception of data from one flight line at 
Guymon (1.6 GHz data on 8/11/78 had highly erratic brightness 
temperatures), brightness temperatures were quite stable. The highly 
variable brightness temperatures Indicated local unmeasured variations 
In the field. Therefore, the following fields at Guymon were deleted 
from further analysis: fields 10, 13, 14, 15 and 17. 

NSOOl/N^S Processing 

Raw data collected on the NASA C-130 were In analog form. The 
data were converted to digital onto 9-track tapes at NASA/Johnson 

Space Center. Included with the surface data were calibration data 
consisting of digital counts from looks at constant radiance targets 
aboard the plane. The calibration data were then used to convert 
digital counts to radiance. 


Ill 









since radiance Is a function of the solar angle, a correction 
factor was needed before comparing crop radiance differences. All the 
Dalhai data were normalized to August 18, the day with the smallest 
solar zenith angle; Guymon data were adjusted to August 11 zenith 
angle conditions. The correction factor used was 

R = — ! (A2) 

cos 6 

where Ri and R^ are the non-normal 1 zed and normalized radiance 
values, respectively, and e Is the solar zenith angle. 

Most of the visible/infrared data were good quality at Dalhart 
and Guymon. One of the exceptions was the excessively noisy water 
absorption bonds (channels 6 and 7) on 8/14/80 at Dalhart. Since no 
means was possible to correct the bands, they were eliminated fro«n 
further data analysis. Also, channel 1 data for fields 6, 8, 10, 12 
and 22 were deleted due to unstable calibration. 

With the exception of channel y (770-863 nm) MMS data at Guymon, 
the calibration Information proved to be quite stable. Table A4a 
lists the equations used to convert raw digital counts to radiance 
values. . Note channel 9 had three different equations applicable at 
different periods of the experiment. 

All of the working NSOOl bands had less stable calibration 
Information at Dalhart. Table A4b lists the equations used to convert 
digital counts to radiance values. Note several channels had 
different calibration values on each flight day. 

Calibration of the thermal band proved to be different for Guy- 
mon and Dalhart. The calibration, using the PRT-5 data, showed that 
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TABLE A4. (Continued) 


ORIGINAL PAGE IS 
OF POOR QUALITY 


5 

5 

5 

6 

6 

6 


„ _ 5.43x10'^* 

^ — 

„ _ 5.43x10’^^ 

^ ~w — 

R * 5.43x10~^^ 

„ 2.8x10"^ * 

^ - T25 


(DC-8) 

(DC-9) 

(DC-9) 

(OC-12) 


(8/14-8/16 (Fit 1)) 
(.3/16 (Fit 2)) 
(8/18) 

(8/16) 


R = ^ *(0C-12) (8/18) 

R = *(DC-16) (8/16 & 8/18) 


at Guymon the low temperature calibration black box aboard the plane 
was too high while the high temperature calibration black box was 
measuring the proper temperature. This implied that low surface 

temperatures were as much as 5 K too high. At Dalhart, the opposite 
condition occurred. The low temperature calibration box was reading 
the proper temperature while the high temperature calibration box was 
reading 5 K too low, suggesting that high surface temperatures were as 
much as 5 K too low. 

The normalization solar correction factors (cos 6i) for Oalhart 
are as follows: August 14, 5.7; August 16 (flight 1), 2.0; and 

(flight 2), 1.1 and August 18, 1.0. For Guymon, the normal ization 

solar correction factors are August 2, 1.7; August 5, 1.6; August 8, 
5.0; August 11, 1.0; August 14, 1.6 and August 17, 1.6. To normalize 
the two data sets, the Guymon data set required a multiplication 

factor of 1.3 to roughly match the radiance values at Oalhart. 

Soil Moisture 

Each sensor has a different cell size. Consequently, to compare 
data, field averages were determined for each sensor and compared to 
ground data field averages. Unfortunately, in some cases, averaging 
point locations of soil moisture proved not to be a reliable field 
average. For instance, several rows were irrigated and seen by the 

sensors but not sampled within the field. Also rainfall events 
occurred at Guymon between sampling periods--on 8/2 and 8/8/78. An 
attempt was made to correct the soil moisture by adding the amount of 
rainfall or irrigation, assuming complete infiltration. In some 
cases, this correction did a good job. But in the end the question- 
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able soil moisture data were deleted from the data set. The fields at 
Guymon with deleted soil moisture data were for 8/2: 22, 27, 20, 25, 
19, 24, 8/8: lx, 2x, 2, 10 and 8/17: lx, (line 2). 

With the deletions, calibrations, and normalizations the Guymon 
and Dalhart data sets were as complete as possible. 
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APPENDIX B 


SCATTERONETER RESPONSES TO SOIL MOISTURE 
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OKlGnMAL PAGE IS 
OF POOR QUALITY 


GUYMON BARE FIELDS 


13.3 W (to DCC) VS VOL. SOR. MOISTURE 


13.3 W (IS DEC) VS VOL. SOL MOISTURE 


? H 




-5- V 


R = 0.43 


R =0.47 


0 S 10 IS 20 2S 30 
VOL. SOIL HOtSTUHC 10-2 CHI 


0 S 10 IS 20 2S 30 
VOL. SOIL noisrunc 10-2 chi 


1 3.3 W (20 OEC) vs VOL. SOL MOISTURE 


13.3 W (40 OEC) VS VOL. SOL MOISTURE 



0 S 10 IS 20 2S 30 
VOL. SOIL HOISruHE 10-2 CHI 



R =0.10 


0 S 10 IS 20 2S 30 
VOL. SOIL HOISruHE 10-2 CHI 


LECENOi riELO x k x CIO 
V V Y C17 
Z Z Z C2 
A 4 A C2I 

• t • CZI4Z6 

• • • C641M 


FIG. B1 Guymon 13.3 GHz VV o° vs. volinnetrlc soil moisture (0*2 cm) 
at 10, 15, 20, and 40 degree Incident angles. 




ORIGINAL PAGE IS 
OF POOR QUALITY 


4.7S HH (10 OCC) VS VOL. SOIL MOBIURC 


GUYMON BARE FIELDS 

4.7S HH (15 OGC) VS VOL. SOIL MOISTURC 



VOL. SOIL noiSTunc io-2 cm 



VOL. SOIL Hoisrunc io-2 chi 


4.75 HH (20 OCC) VS VOL. SOIL M05TURC 


4.75 HH (40 OCC) VS VOL. SOIL MOISTURC 



VOL. SOIL HOISTUflC 10-2 CHI 

UCCNOi MCLO 



VOL. SOIL Noisrunc IO-2 CHI 


X X X CIO 
V Y V 017 
Z Z Z C2 

A 4 A C2I 

• • • C2l4?e 

• • • C641U 


riG. B3 Guymon 4.75 HH o° vs. volumetric soi l"'lfolsture (0-2 cm) at 10, 
15, 20, and 40 degree Incident angles. 
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ORIGINAL PAGE 13 
OF POOR QUALITY 


FIG. 


DALHART BARE FIELDS 


4.7S HN (to OCC) V$ VOL. SO». MOISTURC 



VOL. SOIL nOISTUnC IO-2 CHI 


4.79 HH (15 OCC) VS VOL. SOIL MOISTURC 



VOL. SOIL noiSTunc lo-? cm 


4.75 HH (20 OCC) VS VOL. SOIL MOISTURC 


4.75 HH (40 OCC) VS VOL. SOIL MOISTURC 




VOL. SOIL NOISUmC 10-2 CHI VOL. SOIL HOISTUHI 10-2 CHI 

LCCCNOt FULO • « « 019414 

• * * 0IS4I0 
O O O 017410 

• 4 • 010420 

• • • 021422 
RC6 

B4 Dal hart 4.75 HH o° vs. volumetric soil moisture (0-2 cm) at 
10, 15, 20, and 40 degree Incident angles. 


OPIGJNAL PAQEJS 
OF POOR QUALITY 


GUYMON BARE FIELDS 


4.7A »:V (10 OCC) VS VOL. SOIL MOISTURC 4.7S HV (IS OCC) VS VOL SOIL MOISTURC 



VOL. SOIL MOISTURC 10-2 CHI 



I 

VOL. SOU HOlSTimE \ 0-2 CHI 


4.75 MV (20 DGC) VS VOL. SOIL MOISTURC 



4.75 MV (40 DCC) VS VOL. SOIL MOISTURC 



VOL. SOIL Moisumc C0-? cm 


VOL. SOIL MOISTUMt 10-2 CMI 


LCCCNOi riCLO 


K X K CIO 

V T V 017 
2 Z Z 02 

A A A 021 

• A • C 2 I 426 

• • • 064 m 
MCC 


FIG. B5 Guymon 4.75 HV o” vs. volumetric soil moisture (0-2 cm) at 
10, 15, 20, and 40 degree Incident angles. 





ORIGINAL PAQb 13 
C7 POOR QUALITY 


GUYMON BARE FIELDS 


l.« KH (10 DCC) VS VOL. SOIL MOISTURE 


1.0 HN (IS DEC) VS VOL. SON. MOISTURE 



S -10- 

I : • / 

R 1 




R* -0.48 


R* «n *;ft 


0 s 10 IS 20 25 so 
VOL. SOIL HonTune 10-2 cm 


0 s 10 IS 20 2S so 
VOL. SOIL NOisrune 10-2 cm 


1.0 HH (20 DEC) VS VOL. SOIL MOISTURE 


1.0 HH (40 DEC) VS VOL. SOIL MOISTURE 


8 -12- 
\ 


R* -0.47 


3 -20- * , X 

A- 
• j 


0 S 10 IS 20 25 SO 
VOL. SOIL HOISIU^C 10-2 CMI 


LCceNUi ncLO 


0 5 10 IS 20 2S SO 

VOL. SOIL NOISTU^e 10-2 CHI 


K GIO 
V CI7 
2 C2 

A 021 

• C2U2I 

• 064U 


FIG. B7 Guymon 1.6 GHz HH 0 ° vs. volumetric soil moisture (0-2 cm) 
at 10, 15, 20, and 40 degree Incident angles. 



ORIGINAL PAGE 13 
OF POOR QUALITY 


GUYMON BARE FIELDS 

1.« HV (10 OtC) VS VOl. SOIL MOfSTURC M HV (IS OCC) VS VOl. SOU MOISTUR*’ 



VOL. SOU notsume io-2 cm voi. sou HOisrunc io-2 chi 


1.5 HV (20 OCC) vs VOL. SOIL MOISTURC 1.5 HV (40 OCC) VS VOL. SOIL MOISTURC 



VOl. SOU MOtSTunc 10-2 cm vol. sou HOisrunc io-2 cm 


LCCCNOi riCLO 


H K CIO 
V V 017 
X t. C2 


A A C2I 


A A C2IA2I 


f. • • • C04III 

FIG. B8 Guymon 1.6 GHz HV o” vs. volumetric soil moisture (0-2 cm) 
at 10, 15, 20, and 40 degree Incident angles. 






ORIGINAL FhCE .3 
OF POOR Q 'A'jty 

DALHART BARE FIELDS 


1.« HV (10 DEC) VS VOL. SOIL MOISTUtC 


1.« MV (IS OCC) VS VOL. SOIL MOISTL'RC 


3 - 20 - 

5 



3 -30- ^ 

I 

M 

R 



= 0.02 


R* =0.28 


0 5 10 IS 20 25 30 

VOL. SOIL MOISTURE 10-2 CMI 


0 S 10 IS 20 2S 30 
XOL. SOU MOISTURE 10-2 CNI 


1.1 HV (20 DEC) VS VOL. SOIL MOISTURE 


1.« HV (40 OCC) VS VOL. SOIL MOISTURE 




♦ 




R* =0.36 


R' =0.17 


0 5 10 tS >0 2S SO 

VOL. SOIL MOisrunt lo-? chi 


0 S 10 IS ?0 2S so 
VOL. sou noisrune lo ? cm 


LCCCNOi neic ... ois«i« 

4 ♦ ♦ 0IS«I6 
O O O 017*16 
... 0l6«20 
... 02l«?2 

FIG. B9 Dalhart 1.6 GHz HV o° vs. volumetric soil moisture (0-2 cm) 
at 10, l‘j, 20, and 40 degree Incident angles. 



ORIGINAL PAGE IS 
OF POOR QUALITY 

GUYMON BARE FIELDS 


0.4 HH (10 OCC) VS VOL. SON. MOISTURC 


0 S 10 IS 20 25 90 

VOL. SOIL noiSTunc io-2 cm 

0.4 HH (20 OCC) VS VOL. SOIL MOISTUNC 
-lOH 


0.4 HH (19 OCC) VS VOL. SON. MOISTUNC 


0 5 10 IS 20 2S 90 

VOL. 901L nOISTune 10-2 CHI 

0.4 HH (40 OCC) VS VOL. SON. MOISTUNC 
-25H 


i 

1* 

V 

*4 

1 

2 ’ •• 

• 

R* =0.11 

1 

1 


(d R «0.19 

I ■ ■ " I ' ' I ' 1 1 1 1 1 ...... I .... I 

0 S 10 IS 20 2S 90 
VOL. SOIL HOISruNC i0-2 CHI 

leceNDi riCLO 


H 

0 5 10 IS 20 25 SO 

VOL. SOIL HOlSruAC 10-2 CM) 


n X CIO 
V ▼ 017 
Z Z C2 
A 4 C21 

• • C2I420 

• • C64IH 


FIG. BIO Guymon 0.4 GHz HH o® vs. volumetric soil moisture (0-2 cm) 
at 10. 15, 20. and 40 degree incident angles. 





ORIGINAL PAGE IS 
OF POOR QUALITY 


GUYMON BARE FIELDS 

0.4 HV (10 OCC) VS VOL. SOIL MOISTURE 0.4 HV (IS DEC) VS VOL. SOIL MOISTURE 



0 5 10 IS 20 2S 30 OS 10 IS 20 2S 30 


VOL. SOIL nOlSTUnC 10-2 CNI VOL. SOIL MOISTURE 10-2 CHI 


0.4 HV (20 DEC) VS VOL. SOIL MOISTURE 0.4 HV (40 DEC) VS VOL. SOIL MOISTURE 



0 S 10 IS 20 2S 30 OS 10 IS 20 2S 30 


VOL. SOIL MOISTUMt 10-2 CMI VOL. SOIL MOISTUflt 10-2 CMI 

LECENOt MELO k h « CIO 
V r Y CI7 
Z Z Z C2 
t A A C2I 

• • • C21426 

• • • C641M 

FIG. B12 Guymon 0.4 GHz HV o° vs. volumetric soil moisture (0-2 cm) 
at 10, 15, 20, and 40 degree incident angles. 


ORIGINAL PAGE 13 
OF POOR QUALITY 

DALHART BARE FIELDS 


0.4 HV (10 DEC) VS VOL. SOU UOISTURC 0.4 HV (IS DEC) VS VOL. SOIL MOISTURE 



VOL. SOIL MOISTUW CO-2 CHI VOL. SOIL MOlSTUOt 10-2 CMl 


0.4 HV (20 DEC) VS VOL. SOIL MOISTURE 


0.4 HV (40 DEC) VS VOL. SOIL MOISTURE 



VOL. SOIL HOISTURE 10-2 CHI 

LECCNOi FIELD 



VOL. SOIL HOISTUHE 10-2 CHI 


« * « 01S(I« 

* * * 0IS416 
O O O 017410 

* * » 011420 

* • • D2I422 


FIG, B13 Dalhart 0.4 GHz HV o” vs, volumetric soil moisture (0-2 cm) 


at 10, 15, 20, and 40 degree Incident angles. 
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